A germline chimeric KANK1-DMRT1 transcript derived from a complex structural variant is associated with a congenital heart defect segregating across five generations.

Structural variants (SVs) pose a challenge to detect and interpret, but their study provides novel biological insights and molecular diagnosis underlying rare diseases. The aim of this study was to resolve a 9p24 rearrangement segregating in a family through five generations with a congenital heart defect (congenital pulmonary and aortic valvular stenosis and pulmonary artery stenosis), by applying a combined genomic analysis. The analysis involved multiple techniques, including karyotype, chromosomal microarray analysis (CMA), FISH, genome sequencing (GS), RNA-seq, and optical genome mapping (OGM). A complex 9p24 SV was hinted at by CMA results, showing three interspersed duplicated segments. Combined GS and OGM analyses revealed that the 9p24 duplications constitute a complex SV, on which a set of breakpoints matches the boundaries of the CMA duplicated sequences. The proposed structure for this complex rearrangement implies three duplications associated with an inversion of ~ 2 Mb region on chromosome 9 and a SINE element insertion at the more distal breakpoint. Interestingly, this genomic structure of rearrangement forms a chimeric transcript of the KANK1/DMRT1 loci, which was confirmed by both RNA-seq and Sanger sequencing on blood samples from 9p24 rearrangement carriers. Altogether with breakpoint amplification and FISH analysis, this combined approach allowed a deep characterization of this complex rearrangement. Although the genotype-phenotype correlation remains elusive from the molecular mechanism point of view, this study identified a large genomic rearrangement at 9p24 segregating with a familial congenital heart defect, revealing a genetic biomarker that was successfully applied for embryo selection, changing the reproductive perspective of affected individuals.

Clinical Characterization and Underlying Genetic Findings in Brazilian Patients with Syndromic Microcephaly Associated with Neurodevelopmental Disorders.

Microcephaly is characterized by an occipitofrontal circumference at least two standard deviations below the mean for age and sex. Neurodevelopmental disorders (NDD) are commonly associated with microcephaly, due to perturbations in brain development and functioning. Given the extensive genetic heterogeneity of microcephaly, managing patients is hindered by the broad spectrum of diagnostic possibilities that exist before conducting molecular testing. We investigated the genetic basis of syndromic microcephaly accompanied by NDD in a Brazilian cohort of 45 individuals and characterized associated clinical features, as well as evaluated the effectiveness of whole-exome sequencing (WES) as a diagnostic tool for this condition. Patients previously negative for pathogenic copy number variants underwent WES, which was performed using a trio approach for isolated index cases (n = 31), only the index in isolated cases with parental consanguinity (n = 8) or affected siblings in familial cases (n = 3). Pathogenic/likely pathogenic variants were identified in 19 families (18 genes) with a diagnostic yield of approximately 45%. Nearly 86% of the individuals had global developmental delay/intellectual disability and 51% presented with behavioral disturbances. Additional frequent clinical features included facial dysmorphisms (80%), brain malformations (67%), musculoskeletal (71%) or cardiovascular (47%) defects, and short stature (54%). Our findings unraveled the underlying genetic basis of microcephaly in half of the patients, demonstrating a high diagnostic yield of WES for microcephaly and reinforcing its genetic heterogeneity. We expanded the phenotypic spectrum associated with the condition and identified a potentially novel gene (CCDC17) for congenital microcephaly.

Low-pass whole genome sequencing is a reliable and cost-effective approach for copy number variant analysis in the clinical setting.

INTRODUCTION: Next generation sequencing technology has greatly reduced the cost and time required for sequencing a genome. An approach that is rapidly being adopted as an alternative method for CNV analysis is the low-pass whole genome sequencing (LP-WGS). Here, we evaluated the performance of LP-WGS to detect copy number variants (CNVs) in clinical cytogenetics. MATERIALS AND METHODS: DNA samples with known CNVs detected by chromosomal microarray analyses (CMA) were selected for comparison and used as positive controls; our panel included 44 DNA samples (12 prenatal and 32 postnatal), comprising a total of 55 chromosome imbalances. The selected cases were chosen to provide a wide range of clinically relevant CNVs, the vast majority being associated with intellectual disability or recognizable syndromes. The chromosome imbalances ranged in size from 75 kb to 90.3 Mb, including aneuploidies and two cases of mosaicism. RESULTS: All CNVs were successfully detected by LP-WGS, showing a high level of consistency and robust performance of the sequencing method. Notably, the size of chromosome imbalances detected by CMA and LP-WGS were compatible between the two different platforms, which indicates that the resolution and sensitivity of the LP-WGS approach are at least similar to those provided by CMA. DISCUSSION: Our data show the potential use of LP-WGS to detect CNVs in clinical diagnosis and confirm the method as an alternative for chromosome imbalances detection. The diagnostic effectiveness and feasibility of LP-WGS, in this technical validation study, were evidenced by a clinically representative dataset of CNVs that allowed a systematic assessment of the detection power and the accuracy of the sequencing approach. Further, since the software used in this study is commercially available, the method can easily be tested and implemented in a routine diagnostic setting.

DNA methylation patterns suggest the involvement of DNMT3B and TET1 in osteosarcoma development.

DNA methylation may be involved in the development of osteosarcomas. Osteosarcomas commonly arise during the bone growth and remodeling in puberty, making it plausible to infer the involvement of epigenetic alterations in their development. As a highly studied epigenetic mechanism, we investigated DNA methylation and related genetic variants in 28 primary osteosarcomas aiming to identify deregulated driver alterations. Methylation and genomic data were obtained using the Illumina HM450K beadchips and the TruSight One sequencing panel, respectively. Aberrant DNA methylation was spread throughout the osteosarcomas genomes. We identified 3146 differentially methylated CpGs comparing osteosarcomas and bone tissue samples, with high methylation heterogeneity, global hypomethylation and focal hypermethylation at CpG islands. Differentially methylated regions (DMR) were detected in 585 loci (319 hypomethylated and 266 hypermethylated), mapped to the promoter regions of 350 genes. These DMR genes were enriched for biological processes related to skeletal system morphogenesis, proliferation, inflammatory response, and signal transduction. Both methylation and expression data were validated in independent groups of cases. Six tumor suppressor genes harbored deletions or promoter hypermethylation (DLEC1, GJB2, HIC1, MIR149, PAX6, and WNT5A), and four oncogenes presented gains or hypomethylation (ASPSCR1, NOTCH4, PRDM16, and RUNX3). Our analysis also revealed hypomethylation at 6p22, a region that contains several histone genes. Copy-number changes in DNMT3B (gain) and TET1 (loss), as well as overexpression of DNMT3B in osteosarcomas provide a possible explanation for the observed phenotype of CpG island hypermethylation. While the detected open-sea hypomethylation likely contributes to the well-known osteosarcoma genomic instability, enriched CpG island hypermethylation suggests an underlying mechanism possibly driven by overexpression of DNMT3B likely resulting in silencing of tumor suppressors and DNA repair genes.

Establishment of iPSC lines and zebrafish with loss-of-function AHDC1 variants: Models for Xia-Gibbs syndrome.

Xia-Gibbs syndrome (XGS) is a syndromic form of intellectual disability caused by heterozygous AHDC1 variants, but the pathophysiological mechanisms underlying this syndrome are still unclear. In this manuscript, we describe the development of two different functional models: three induced pluripotent stem cell (iPSC) lines with different loss-of-function (LoF) AHDC1 variants, derived by reprogramming peripheral blood mononuclear cells from XGS patients, and a zebrafish strain with a LoF variant in the ortholog gene (ahdc1) obtained through CRISPR/Cas9-mediated editing. The three iPSC lines showed expression of pluripotency factors (SOX2, SSEA-4, OCT3/4, and NANOG). To verify the capacity of iPSC to differentiate into the three germ layers, we obtained embryoid bodies (EBs), induced their differentiation, and confirmed the mRNA expression of ectodermal, mesodermal, and endodermal markers using the TaqMan hPSC Scorecard. The iPSC lines were also approved for the following quality tests: chromosomal microarray analysis (CMA), mycoplasma testing, and short tandem repeat (STR) DNA profiling. The zebrafish model has an insertion of four base pairs in the ahdc1 gene, is fertile, and breeding between heterozygous and wild-type (WT) animals generated offspring in a genotypic proportion in agreement with Mendelian law. The established iPSC and zebrafish lines were deposited on the hpscreg.eu and zfin.org platforms, respectively. These biological models are the first for XGS and will be used in future studies that investigate the pathophysiology of this syndrome, unraveling its underlying molecular mechanisms.

Waardenburg Syndrome: The Contribution of Next-Generation Sequencing to the Identification of Novel Causative Variants.

Waardenburg syndrome (WS) is characterized by hearing loss and pigmentary abnormalities of the eyes, hair, and skin. The condition is genetically heterogeneous, and is classified into four clinical types differentiated by the presence of dystopia canthorum in type 1 and its absence in type 2. Additionally, limb musculoskeletal abnormalities and Hirschsprung disease differentiate types 3 and 4, respectively. Genes PAX3, MITF, SOX10, KITLG, EDNRB, and EDN3 are already known to be associated with WS. In WS, a certain degree of molecularly undetected patients remains, especially in type 2. This study aims to pinpoint causative variants using different NGS approaches in a cohort of 26 Brazilian probands with possible/probable diagnosis of WS1 (8) or WS2 (18). DNA from the patients was first analyzed by exome sequencing. Seven of these families were submitted to trio analysis. For inconclusive cases, we applied a targeted NGS panel targeting WS/neurocristopathies genes. Causative variants were detected in 20 of the 26 probands analyzed, these being five in PAX3, eight in MITF, two in SOX10, four in EDNRB, and one in ACTG1 (type 2 Baraitser-Winter syndrome, BWS2). In conclusion, in our cohort of patients, the detection rate of the causative variant was 77%, confirming the superior detection power of NGS in genetically heterogeneous diseases.

Burden of Rare Copy Number Variants in Microcephaly: A Brazilian Cohort of 185 Microcephalic Patients and Review of the Literature.

Microcephaly presents heterogeneous genetic etiology linked to several neurodevelopmental disorders (NDD). Copy number variants (CNVs) are a causal mechanism of microcephaly whose investigation is a crucial step for unraveling its molecular basis. Our purpose was to investigate the burden of rare CNVs in microcephalic individuals and to review genes and CNV syndromes associated with microcephaly. We performed chromosomal microarray analysis (CMA) in 185 Brazilian patients with microcephaly and evaluated microcephalic patients carrying < 200 kb CNVs documented in the DECIPHER database. Additionally, we reviewed known genes and CNV syndromes causally linked to microcephaly through the PubMed, OMIM, DECIPHER, and ClinGen databases. Rare clinically relevant CNVs were detected in 39 out of the 185 Brazilian patients investigated by CMA (21%). In 31 among the 60 DECIPHER patients carrying < 200 kb CNVs, at least one known microcephaly gene was observed. Overall, four gene sets implicated in microcephaly were disclosed: known microcephaly genes; genes with supporting evidence of association with microcephaly; known macrocephaly genes; and novel candidates, including OTUD7A, BBC3, CNTN6, and NAA15. In the review, we compiled 957 known microcephaly genes and 58 genomic CNV loci, comprising 13 duplications and 50 deletions, which have already been associated with clinical findings including microcephaly. We reviewed genes and CNV syndromes previously associated with microcephaly, reinforced the high CMA diagnostic yield for this condition, pinpointed novel candidate loci linked to microcephaly deserving further evaluation, and provided a useful resource for future research on the field of neurodevelopment.

Chromosomal microarray analyses from 5778 patients with neurodevelopmental disorders and congenital anomalies in Brazil.

Chromosomal microarray analysis (CMA) has been recommended and practiced routinely since 2010 both in the USA and Europe as the first-tier cytogenetic test for patients with unexplained neurodevelopmental delay/intellectual disability, autism spectrum disorders, and/or multiple congenital anomalies. However, in Brazil, the use of CMA is still limited, due to its high cost and complexity in integrating the results from both the private and public health systems. Although Brazil has one of the world's largest single-payer public healthcare systems, nearly all patients referred for CMA come from the private sector, resulting in only a small number of CMA studies in Brazilian cohorts. To date, this study is by far the largest Brazilian cohort (n = 5788) studied by CMA and is derived from a joint collaboration formed by the University of São Paulo and three private genetic diagnostic centers to investigate the genetic bases of neurodevelopmental disorders and congenital abnormalities. We identified 2,279 clinically relevant CNVs in 1886 patients, not including the 26 cases of UPD found. Among detected CNVs, the corresponding frequency of each category was 55.6% Pathogenic, 4.4% Likely Pathogenic and 40% VUS. The diagnostic yield, by taking into account Pathogenic, Likely Pathogenic and UPDs, was 19.7%. Since the rational for the classification is mostly based on Mendelian or highly penetrant variants, it was not surprising that a second event was detected in 26% of those cases of predisposition syndromes. Although it is common practice to investigate the inheritance of VUS in most laboratories around the world to determine the inheritance of the variant, our results indicate an extremely low cost-benefit of this approach, and strongly suggest that in cases of a limited budget, investigation of the parents of VUS carriers using CMA should not be prioritized.

Genetic investigation of syndromic forms of obesity.

BACKGROUND: Syndromic obesity (SO) refers to obesity with additional phenotypes, including intellectual disability (ID)/developmental delay (DD), dysmorphic features, or organ-specific abnormalities. SO is rare, has high phenotypic variability, and frequently follows a monogenic pattern of inheritance. However, the genetic etiology of most cases of SO has not been elucidated. SUBJECTS AND METHODS: In this study, we investigated 20 SO patients by whole-exome sequencing (WES) trios to identify causal genetic variants. RESULTS: 4/20 patients had negative results for array comparative genomic hybridization (aCGH) analyses. In the remaining 15 patients, in addition to SNVs and indels, CNVs were also evaluated. Pathogenic/likely pathogenic (P/LP) SNVs/indels were detected in 6/20 patients (involving MED13L, AHDC1, EHMT1, MYT1L, GRIA3, and SETD1A), while two patients carried an inherited VUS. In addition, P/LP CNVs were observed in 3/15 patients (involving SATG2, KIAA0442, and MEIS2). CONCLUSIONS: All nine detected P/LP variants involved genes already known to lead to syndromic ID/DD; however, for only two genes (EHMT1 and MYT1L) is the link with obesity well established. This is the first study applying a comprehensive genomic investigation of an SO cohort, showing a high diagnostic yield (~47%). Additionally, our findings suggested that several known ID/DD genes may also predispose individuals to SO.

Copy number variations in a Brazilian cohort with autism spectrum disorders highlight the contribution of cell adhesion genes.

Prediction of pathogenicity of rare copy number variations (CNVs), a genomic alteration known to contribute to the etiology of autism spectrum disorder (ASD), represents a serious limitation to interpreting genetic tests, particularly for genetic counseling purposes. Chromosomal microarray analysis (CMA) was conducted in a unique collection of 144 Brazilian individuals with ASD of strong European and African ancestries. Rare CNVs were detected in 39 patients: 41 of unknown significance (VUS), four pathogenic and one likely pathogenic CNVs (clinical yield of 4.1%; 5/122). Based on gene content and recurrence in three large cohorts [a Brazilian neurodevelopmental disorder cohort, the autism MSSNG cohort, and the Canadian-based Centre for Applied Genomics microarray database], this work strengthened the pathogenicity of 14 genes (FAT1, CAMK4, BIRC6, DPP6, CSMD1, CTNNA3, CDH8/CDH11, CDH13, OR1C1, CNTN6, CNTNAP4, FGF2 and PTPRN2) within 14 CNVs. Notably, enrichment of cell adhesion proteins to ASD etiology was identified (p < 0.05), highlighting the importance of these gene families in the etiology of ASD.

A Small Supernumerary Xp Marker Chromosome Including Genes NR0B1 and MAGEB Causing Partial Gonadal Dysgenesis and Gonadoblastoma.

Copy number variations of several genes involved in the process of gonadal determination have been identified as a cause of 46,XY differences of sex development. We report a non-syndromic 14-year-old female patient who was referred with primary amenorrhea, absence of breast development, and atypical genitalia. Her karyotype was 47,XY,+mar/46,XY, and FISH analysis revealed the X chromosome origin of the marker chromosome. Array-CGH data identified a pathogenic 2.0-Mb gain of an Xp21.2 segment containing NR0B1/DAX1 and a 1.9-Mb variant of unknown significance from the Xp11.21p11.1 region. This is the first report of a chromosomal microarray analysis to reveal the genetic content of a small supernumerary marker chromosome detected in a 47,XY,+der(X)/46,XY karyotype in a non-syndromic girl with partial gonadal dysgenesis and gonadoblastoma. Our findings indicate that the mosaic presence of the small supernumerary Xp marker, encompassing the NR0B1/DAX1 gene, may have been the main cause of dysgenetic testes development, although the role of MAGEB and other genes mapped to the Xp21 segment could not be completely ruled out.

Congenital chromoanagenesis in the routine postnatal chromosomal microarray analyses.

Chromosomal microarray analyses (CMA) have greatly increased both the yield and diagnostic accuracy of postnatal analysis; it has been used as a first-tier cytogenetic test in patients with intellectual disability, autism spectrum disorder, and multiple congenital abnormalities. During the last 15 years, we performed CMA in approximately 8,000 patients with neurodevelopmental and/or congenital disorders, of which 13 (0.16%) genetically catastrophic complex chromosomal rearrangements were identified. These ultrarare rearrangements showed clustering of breakpoints, characteristic of chromoanagenesis events. Al1 13 complex events display underlying formation mechanisms, originating either by a synchronization of the shattering of clustered chromosome regions in which regional asynchrony of DNA replication may be one of the main causes of disruption. We provide an overview of the copy number profiling in these patients. Although several previous studies have suggested that chromoanagenesis is often a genetic disease source in postnatal diagnostic screening, due to either the challenge of clinical interpretation of these complex rearrangements or the limitation of microarray resolution relative to the small size and complexity of chromogenic induced chromosome abnormalities, bringing further attention and to study its occurrence in the clinical setting is extremely important.

Detection of mosaicism for segmental and whole chromosome imbalances by targeted sequencing.

Mosaic segmental and whole chromosome copy number alterations are postzygotic variations known to be associated with several disorders. We have previously presented an efficient targeted sequencing approach to simultaneously detect point mutations and copy number variations (CNVs). In this study, we evaluated the efficiency of this approach to detect mosaic CNVs, using seven postnatal and 19 tumor samples, previously characterized by chromosomal microarray analyses (CMA). These samples harbored a total of 28 genomic imbalances ranging in size from 0.68 to 171 Mb, and present in 10-80% of the cells. All CNV regions covered by the platform were correctly identified in postnatal samples, and only seven out of 19 CNVs from tumor samples were not identified either because of a lack of target probes in the affected genomic regions or an absence of minimum reads for an alteration call. These results demonstrate that, in a research setting, this is a robust approach for detecting mosaicism in cases of segmental and whole chromosome alterations. Although the current sequencing platform presented a resolution similar to genomic microarrays, it is still necessary to further validate this approach in a clinical setting in order to replace CMA and sequencing analyses by a single test.

Manifesting carriers of X-linked myotubular myopathy: Genetic modifiers modulating the phenotype.

OBJECTIVE: To analyze the modulation of the phenotype in manifesting carriers of recessive X-linked myotubular myopathy (XLMTM), searching for possible genetic modifiers. METHODS: Twelve Brazilian families with XLMTM were molecularly and clinically evaluated. In 2 families, 4 of 6 and 2 of 5 manifesting female carriers were identified. These females were studied for X chromosome inactivation. In addition, whole-exome sequencing was performed, looking for possible modifier variants. We also determined the penetrance rate among carriers of the mutations responsible for the condition. RESULTS: Mutations in the MTM1 gene were identified in all index patients from the 12 families, being 4 of them novel. In the heterozygotes, X chromosome inactivation was random in 3 of 4 informative manifesting carriers. The disease penetrance rate was estimated to be 30%, compatible with incomplete penetrance. Exome comparative analyses identified variants within a segment of 4.2 Mb on chromosome 19, containing the killer cell immunoglobulin-like receptor cluster of genes that were present in all nonmanifesting carriers and absent in all manifesting carriers. We hypothesized that these killer cell immunoglobulin-like receptor variants may modulate the phenotype, acting as a protective factor in the nonmanifesting carriers. CONCLUSIONS: Affected XLMTM female carriers have been described with a surprisingly high frequency for a recessive X-linked disease, raising the question about the pattern of inheritance or the role of modifier factors acting on the disease phenotype. We demonstrated the possible existence of genetic mechanisms and variants accountable for the clinical manifestation in these women, which can become future targets for therapies.

Investigating Genetic Factors Contributing to Variable Expressivity of Class I 17p13.3 Microduplication.

17p13.3 microduplications are rare copy number variations (CNVs) associated with variable phenotypes, including facial dysmorphism, developmental delay, intellectual disability, and autism. Typically, when a recognized pathogenic CNV is identified, other genetic factors are not considered. We investigated via whole-exome sequencing the presence of additional variants in four carriers of class I 17p13.3 microduplications. A 730 kb 17p13.3 microduplication was identified in two half-brothers with intellectual disability, but not in a third affected half-brother or blood cells from their normal mother (Family A), thus leading to the hypothesis of maternal germline mosaicism. No additional pathogenic variants were detected in Family A. Two affected siblings carried maternally inherited 450 kb 17p13.3 microduplication (Family B); the three carriers of the microduplication exhibited microcephaly and learning disability/speech impairment of variable degrees. Exome analysis revealed a variant of uncertain significance in RORA, a gene already linked to autism, in the autistic boy; his sister was heterozygous for a CYP1B1 pathogenic variant that could be related to her congenital glaucoma. Besides, both siblings carried a loss-of-function variant in DIP2B, a candidate gene for intellectual disability, which was inherited from their father, who also exhibited learning disability in childhood. In conclusion, additional pathogenic variants were revealed in two affected carriers of class I 17p13.3 microduplication (Family B), probably adding to their phenotypes. These results provided new evidence regarding the contribution of RORA and DIP2B to neurocognitive deficits, and highlighted the importance of full genetic investigation in carriers of CNV syndromes with variable expressivity. Finally, we suggest that microcephaly may be a rare clinical feature also related to the presence of the class I 17p13.3 microduplication.

Mechanistic insights revealed by a UBE2A mutation linked to intellectual disability.

Ubiquitin-conjugating enzymes (E2) enable protein ubiquitination by conjugating ubiquitin to their catalytic cysteine for subsequent transfer to a target lysine side chain. Deprotonation of the incoming lysine enables its nucleophilicity, but determinants of lysine activation remain poorly understood. We report a novel pathogenic mutation in the E2 UBE2A, identified in two brothers with mild intellectual disability. The pathogenic Q93E mutation yields UBE2A with impaired aminolysis activity but no loss of the ability to be conjugated with ubiquitin. Importantly, the low intrinsic reactivity of UBE2A Q93E was not overcome by a cognate ubiquitin E3 ligase, RAD18, with the UBE2A target PCNA. However, UBE2A Q93E was reactive at high pH or with a low-pKa amine as the nucleophile, thus providing the first evidence of reversion of a defective UBE2A mutation. We propose that Q93E substitution perturbs the UBE2A catalytic microenvironment essential for lysine deprotonation during ubiquitin transfer, thus generating an enzyme that is disabled but not dead.